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Abstract 

The general expression of the Stern-Gerlach force is deduced for a relativistic 
charged spin-^ particle which travels inside a time varying magnetic field. This 
result was obtained either by means of two Lorentz boosts or starting from 
Dirac's equation. Then, the utilization of this interaction for attaining the spin 
states separation is reconsidered in a new example using a new radio-frequency 
arrangement. 

1 The Relativistic Stern-Gerlach Force 

The time varying Stern-Gerlach, SG, interaction of a relativistic fermion with an 
e.m. wave has been proposed to separate beams of particles with opposite spin states 
corresponding to different energies [1]. We will show how spin polarized particle will 
exchange energy with the electromagnetic field of an RF resonator. 

Let us denote with (x, y, z) the coordinates of a particle in the laboratory, and 
with (x',y',z') the coordinates in the particle rest frame, PRF. In the latter the SG 
force that represents the action of an inhomogeneous magnetic field on a particle 
endowed with a magnetic moment ft is 

fsa = V V " B') = £r(P* ■ B')x + £r(P*- B')y + JL{j? • B')z (1) 

with 

P = 9^S. (2) 

Here e is the elementary charge with " + " for protons and positrons, p, e + , and " =" 
for antiprotons and electrons, p, e~ , making pi and S either parallel or antiparallel to 
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each other, respectively, m is the rest mass of the particle, g the gyromagnetic ratio 
and a the anomaly defined as 



.9-2 
2 



1.793 (0 = 5.586) for p,p 
1.160 x 1(T 3 for e ± 



(3) 



Notice that in Eq.(TjQ) we have defined the magnetic moment as fx* in the rest 
frame, rather than as fjf. In the rest frame the quantum vector S, or spin, has 



modulus 15*1 = ys(s + l)Ti and its component parallel to the magnetic field lines can 
only take the following values 



s, 



s + 1, s — 1, s)h, 



where Ti is the reduced Planck's constant. Combining Eqs.© and 
the magnetic moment in the PRF 

1.41 x 1(T 26 JT- 1 
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9.28 x 10" 24 JT 
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we obtain for 
(5) 
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For a particle traveling along the axis z, the Lorentz transformations of the dif- 
ferential operators and of the force yield 



_d_ = d_ _d_ _ d_ _d_ 

dx 1 dx dy' dy dz 

f± = h f± f\\ — f\\ (fz = f'z) 



' I \dz c dt) 



The force (CQ) is boosted to the laboratory system as 

fsa = X ~UF ■ + • B')y + £(/T • B')z- 

Because of the Lorentz transformation of the fields [3] E, B and E', B' 

E' = 1 (E + cPxB)-^P0-E) 

B' = 1 ( K B-ixEj-^ I 00.B) 

the energy in the rest frame (/7* • B') becomes 

• B') = 7/ i* (B x + &E V ) + j»* y (B y - + 



(6) 
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Combining Eqs.flH} and ([7]), by virtue of Eq.(jSJ), after some algebra we can finally 
obtain the SG force components in the laboratory frame: 



( f - „* f d- 5 * _l §.dEy_\ , * /aBy _ PdEx\ , I *dB z 
•> x ^x \ dx ^~ c dx J ^~ ry \ dx c dx J 7^2 ax 

/ = ^ (ajk + + /i! f ^ - ^ + 

r~x \ dy c dy J <y \ dy c ay J ' jr~z dy 

, /Z — fJ^C ZX + H*yC Z y "I" /i^C^, 
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with 



r a 
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( dB x 
\ dz 



+ 



PdB, 



dt 



+ f( 



BEy (3 dEy 

dz ^ c dt 



a 



Y 



( 9 B y §_8By\ _ §_ ( dE x , P dE x 
\ dz * c dt ) c \ dz ' c at 
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These results can also be obtained from the quantum relativistic theory of the 
spin- 1 charged particle [2]. Let us introduce the Dirac Hamiltonian 



H = ecp + ca ■ {p — eA) + •y mc 2 
having made use of the Dirac 's matrices 
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O a 
-a O 



7o 



J O 
—J 



a = 7o7 



O a 
a O 



where a is a vector whose components are the Pauli's matrices 



-i 

1 



1 
-1 



1 

1 
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I is the 2x2 identity matrix, O the null matrix and having chosen the y-axis parallel to 
the main magnetic field. A standard derivation leads to the non relativistic expression 
of the Hamiltonian exhibiting the SG interaction with the "normal" magnetic moment 



H 



(15) 



which coincides with the Pauli equation and is valid in the PRF. 

To complete the derivation we must add the contribution from the anomalous 
magnetic moment to the SG energy term in the previous equation, with a factor 
1 + a = |, yielding 



(16) 



In order to obtain the z-component of the SG force in the Laboratory frame along 
the direction of motion of the particle, we must boost the whole Pauli term of Eq. (fl5T) 
by using the unitary operator U in the Hilbert space [I], which expresses the Lorentz 
transformation 

U' 1 [s£(7o<? ■ B')] U = g^( 7o a ■ B>) [S-\ lo a x )S + S'\ lo a y )S + S-\ lo a z )S] 

(17) 

that can be written in terms of the equivalent transformation in the 4x4 spinor space 



S = exp 1 70(7 •■£)!} = cosh I + 



a z 
a?. 



sinh I 



with 

v = 



coshw 



7 = Lorentz factor, sinh u = a/7 2 — 1 



(18) 



7) ■ (I 9 ) 
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From Eqs. flIT]) and (|T8|) . due to the algebraic structure of the 7 and a matrices, 
we obtain in the laboratory frame the three components of the SG force 



with 



S -1 (700-3.) S = J (T X 

S^do^S = j Q a y , (20) 

S- 1 (~f <?z)S = 7(700-*) + Ho^Vi 2 - 1 



O X 

75 = IxlylzlO = t[ z Q ) ■ (21) 



From Eqs.( l20|) we can deduce the expectation values of the SG force in the Lab- 
oratory system with a defined spin -along the y-axis in our case- via the expectation 
values of the Pauli matrices and of the Pauli interaction term of the proper force 

f, = wV [(* + If 1 )] ■ (22) 

In our case only the second of Eqs. (!20|) gives a non vanishing result, while both the 
first and third produce a null contribution to the force, because of the orthogonality 
of the two spin states s = ±| and the properties of the a matrices. 



2 The radio-frequency system 

Let us consider the standing waves built up inside a rectangular radio-frequency 
resonator, tuned to a generic TE Modefl]. Resonator dimensions are: width a, height 
b and length d, as shown in FigHJ On the cavity axis, which coincides with the beam 



axis, the electric and 


magnetic fields are [5] 


E x 


= -So(T)#cos(^)sin(^)sin(^)sin^ 


Ey 


= So(^)^sin(^)co S (^)sin(^)sin^ 


E z 


= [as typical for any TE mode] 


B x 


= -it m (?) sin m c ° s (T) c ° s (*¥■) c ° s 


By 


= -If (T) (S) «» m sin (=pt) cos cos art 


B z 


= B oC os(H^)cos(^p)sin(^)cos^ 


where Bq is the RF peak magnetic field, m, n and p are integer mo 



K c = ^) 2 +(f) 2 - (23) 
The angular frequency of the e.m. wave from the RF generator is 



LO 



^RF = ft = <V M + W + (f ) • (24) 
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Figure 1: Sketch of the rectangular cavity. The coordinates of the beam axis are 
x = | and y = |. 



In contrast with an open waveguide, in a bounded cavity we can define a phase 
velocity Vph and a cavity wavelength A wg , as typical of any e.m. in a refractive media, 
according to the relations 



; ph 

c 



and 

It is also 



A w g — Pph^RF- (26) 

V ph = (3 ph c = (3 ph ^ = ^ (27) 
trf trf 

Notice that /3 p h can take any value, even larger than one, since it is freely depen- 
dent on the cavity geometrical parameters. Moreover, combining Eqs. (1241) and ff25l) 
we obtain 

d = |p/5phA RF = |pA wg (28) 



which describes the connection between the cavity length d and the wavelengths, as 
shown in Figj2j For simplicity, let's choose the transverse electric mode TEoi p , so 
Eqs. ( 124]) and (T25]) reduce respectively to 

^ = ^RF = ^(f) 2 + (f) 2 = /3 P < and fa = sjl+(ff ( 29 ) 
or, setting the mode index p = 1, 

^ = ^RF = cy(f) 2 + (i) 2 = &< and /? ph = yi + (f) 2 , (30) 
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Figure 2: Vertical component (inside the cavity) and fringes (at both cavity ends) of 
By for p = 1. 

which are the quantities pertaining to the preferred TE n mode whose non zero field 
components on the cavity axis are 

f B y (z,t) = --Bo|cos(^)cosu;t 

(31) 

[ E x (z, t) = -uB ± sin (f ) smut 

It is important to emphasize that in all the field components met so far there is 
a clear separation between spatial and temporal contributions, as typical of standing 
waves. Besides, the boundary conditions of the electric and magnetic fields of the 
e.m. dictate the shape of the spatial component which, in turn, oscillates in time 
with the frequency ujrf. Then, at the cavity entrance and exit the field components 
( I3T]) become on axis 

( B y (0,t) = -Bo^coscot 
Entrance =>- < . (32) 

{ E x (0,t) =0 

and 

!B y (d, t) = —B 2 cos7tcoscg>£ = £>o^ cosu;£ 
(33) 
E x (d, t) = —u^ sin7rsinu;t = 

where t is a generic time. The null values of E x at the cavity ends confirm a typical 
pattern of the transverse electric mode. 
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Figure 3: Edge fields at both ends of a single cavity for p equal to either an even 
number (left) or an odd number (right). 

3 Stern-Gerlach interaction with the cavity field 

From Eq. (l22p . after some algebra, we obtain tha a charged fermion which crosses a 
radio-frequency resonator, tuned on the TE n mode, acquires (or loses) an energy 
amount when interacts with the field component in the "body" of the cavity shown 
in Fig. EC] 

(^Ox-tag = J q f z dz = J q „ C zy dz = /?7% ~ d ^ _ p2 (I + cos — ttJ 

(34) 

still assuming that the spin is not precessing. 

However, since the cavity cannot be completely enclosed but must have apertures 
at both ends to allow the particle bean to pass through and consequently will have 
fringe fields, in order to calculate the full SG interaction it is necessary to deal with 
the interaction with these fields. This is discussed right below. 

3.1 Fringe fields 

In order to fulfill the boundary conditions f[3"2"j) and ([3"3"j) . a cavity tuned in its TE 011 
mode must be exactly filled by either an even or an odd number of cavity dependent 
half wave-lengths, Eq. (l26|) . as illustrated in Figs. [2] and [31 

Consider now a bunch of particles crossing the cavity in synchronism with the RF 
field. This requires that the bunch centre of mass that enters the cavity at the instant 
t = and would leave the cavity at t = trf, at magnetic field values, respectively 

5 tf (0,0) = -B | and B y (d, r RF ) = (35) 

The field values at both ends fade rapidly to zero over a small distance \8\ just 
outside the cavity (see Figures.) We may consider these fringe fields as small- valued 
functions in the (y, z)-plane, since the time St necessary for a particle to proceed 
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through this distances can be very small in comparison with trf, depending of course 
by the size of the beam channel, or 

i B vk = -Mai*) with g(-6)=0, g(0) = 1 

(36) 

[A/Lt = B oih(z) with h(d) = l, h(d + 6) = 

Under these conditions, a relativistic fermion with its spin directed along the y-axis 
and traversing the cavity will experience a SG force parallel to the z-axis (direction 
of motion), see Eq. tfTOl 

/* = fCty (37) 

where C zy is given by the second of the set of Eqs.f llll) . For the moment we assume 
that the spin will conserve its orientation during traversal 

The electric field E x and its derivatives in this equation are almost constantly 
zero, because of the boundary conditions on the walls of the cavity and at the extreme 
points z = —5 and z — d + 8. Furthermore, the function^^ 2 -^ is almost zero along 
the fringe segments because of its proportionality to smut, with t equal to the 5t 
mentioned before. Consequently we have 

C zy ^l 29 -t, (38) 

and for the entire fringe field 

(39) 

{ t/Jout = (^) 

Making use of eqs. (1351) and f[3"9~j) . the energy increments [SU]- m and [5£/] ou t related 
to the fringe fields are easily evaluated since the integrals J® 5 f z dz and f$ +s f z dz only 
depend upon the extreme points (|36|) and do not depend on the curve that connects 
them. In fact f z dz becomes an exact differential. Then we obtain for the energy 
exchange at both edges 

(5U) in = (5U) out = -B fi*^ 2 . (40) 
The total energy exchange at the edges is therefore 

(SU)ft = (^)in + («0aut = -2B 0f I*^ 2 . (41) 

3.2 Full energy interaction 

By adding the fringe contributions (j4"Ij) to the cavity body crossing contribution (jMj) 
seen before, obtain 

(SU) tot = (5U) S + (5U) X = - 7 2 5oM^/(/3 ph , (3) (42) 
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with 



9 02 ^(2-^)-! 



(l + COS ^7T 

For ultra relativistic particles (j3 ~ 1) Eq. (142]) reduces to 

(5Z7)tot ^ -7 2 5 /i*^(l - cos/9 ph 7r). 
This last result deserves a few comments. In fact, if we set 



P h 



±A 



2 W S 



A 



RF 



(43) 



(44) 



(45) 



the total energy contribution ( T441) vanishes, implying a full cancellation of the effect. 
On the other hand if we set 



(3 P h — 3 

the total energy contribution (1441) becomes 

(5[/)tot ^ -2 7 2 5 /x 



d — |A wg — |Arf 



(46) 



(47) 



as deduced from Eq. fl28|) . In Table I we gather values calculated from Eq. (j43~]) for 
non-relativistic and ultra-relativistic particles for, either j3 p h = 2 or /3 p h = 3 at two 
proton energies. Each /3 p h is accompanied by the corresponding ratio cavity-length 
over cavity-height. 

Table I: f(f3 ph ,/3) 



(3 P h => | 


Low Energy 
(e.g. W kiD = 5 MeV) 


High Energy 
(e.g. W kin = 30 GeV) 


2 => 1.732 


2.01 





3 => 2.828 


2.02 


2 



Furthermore, if we consider two contiguous cavities, there will be a gradient be- 
tween the positive B y at the end of the first cavity and a negative B y at the beginning 
of the second cavity, as shown in Fig. |H In this case we may consider the magnetic 
field at the interface as linearly dependent on z, that is 



Reiterating what done before, obtain 



-2B ±z. 



(4* 



/, = -2SoM*^7 2 (49) 
(<5(7)cav2cav = (SU) CC = /j, f,dz = -2B 0j u*^ J[| - (- 1)] ^ = -2£ M* \l 2 
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z 



d 



n-th cavity (n+l)-th cavity 

Figure 4: Magnetic field gradient between two contiguous cavities 

which means that, for N cavities, we shall have as final result for ultra relativistic 
particles 

r o f 0r/ v = 2 

{8U) tot = N(6U) X -in e -(N-l)(6U) cc -(6U) g = l . (50) 

I li B oM*7 2 for& h = 3 

Conversely, if /3 p h is even, particles with their spin pointing always in the same 
direction cannot exchange energy with the standing wave of a TE resonator. A spin 
rotator j6] can align the particle magnetic moments either parallel or anti-parallel 
to the directions of the magnetic field gradients, thus allowing the desired energy 
interaction. This situation would be similar to what happens in a multi-stage tandem 
van de Graaff, where the ions are repetitively accelerated by the same electrostatic 
field, becoming alternatively negative, via an addition of electrons, or positive, via 
electron stripping. 

Unfortunately, the field integral (B^d = fin 1 ^- = (3 5.46 Tm, for p, p) for attaining 
a spin rotation is so large that this solution is unpractical. Instead, the example 
of /? p h equal to an odd number seems much more suitable since does not require 
cumbersome magnets, but only longer cavities (compare Eqs. (j4"5l) and (j46l) ). In 
fact, the magnetic moments are (de) accelerated by the field tails at the cavity ends, 
while don't change their energy when crossing the cavities. This situation resembles 
the Wideroe linac where the charged particles are accelerated by the electric fields 
between two contiguous drift tubes, but don't change their energy while crossing the 
tubes themselves. 
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4 Concluding remarks 



On the basis of the previous estimates, we feel ready to propose the time varying SG 
interaction as a method for attaining a spin state separation of an unpolarized beam 
of, say (anti)protons, since the energy of particles with opposite spin orientations will 
differ and beams in the two states can be separated. In a first stage of the study of a 
sensible practical design, we intend to proceed with numerical simulations. As a first 
step, we intend to verify the correctness of Eqs.( j42l and (pf3l) setting once /? p h = 2 and 
then /? p h = 3, in a cavity where the field line pattern can be realistically controlled. 

Beyond the verification of the present theory, there is also the aim of studying the 
effects generated by the spin precession inside the cavity, that we did not yet address 
in this note. 

Next, we shall consider a spin splitter scheme based on the lattice of an existing or 
planned (anti)proton ring endowed with an array of splitting cavities. The principal 
aim of the latter implementations is to check the mixing effect |7]|8] of the longitudinal 
phase-plane filament at ion, i.e. the actual foe which could frustrate the entire spin 
splitting process. 
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